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We use a Landau theory and the e-expansion to study the superfluid phase transition of two-color 
QCD at nonzero temperature, T, and baryonic chemical potential, fj,. At low T, and for Nf flavors of 
massless quarks, the global SU(N/) xSU(N/) xU(l) symmetry is spontaneously broken by a diquark 
condensate down to Sp(N/) xSp(N/) for any fi > 0. As the temperature increases, the diquark 
■ condensate melts, and at sufficiently large T the symmetry is restored. Using renormalization group 

' arguments, we find that in the presence of the chiral anomaly term there can be a second order 

phase transition when N/ = 2 or N/ > 6, while the transition is first order for N/ = 4. We discuss 
the relevance of these results for the emergence of a tricritical point recently observed in lattice 
simulations. 
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I. INTRODUCTION 

(N , 

During the last few years, a considerable amount of interest has been devoted to the phases of QCD at a nonzero 
baryonic chemical potential /i. Based on either one-gluon exchange or instanton induced interactions, there are 
compelling theoretical indications that a direct attractive quark-quark channel gives rise to diquark condensation at 
asymptotically large fi, similar in spirit to the BCS-scenario [1]. In particular, very dense quark matter is unstable 
y—i \ to the formation of a diquark condensate, and thus is generally believed to be a color superconductor at large /z. 
P\| 1 The QCD-based predictions rely on perturbation theory and are trustworthy when the coupling constant is small, 
g([i) <C 1, i.e. (i » Aqcd, where Aqcd is the hadronic scale (see [2]). To analytically study the properties of 
QCD at more moderate densities, where the quark chemical potential is only a few times Aqcd, one has to rely on 
fS ,' other methods. For instance, one can resort to phenomenological models that are inspired by, but not derived from, 
QCD. An alternative approach would be to use numerical simulations, since lattice gauge theory is a first-principle 
method capable of making quantitative predictions of strong-coupling QCD. At T > and /i = 0, the numerical 
results show that there is a phase transition between the confined and deconfined phases of pure glue matter, at a 
critical temperature T c ~ 270 MeV [3]. With massless quarks, there is a transition to a chirally symmetric phase at 
T c ~ 155 — 175 MeV [4], and recent simulations also show that deconfinement occurs simultaneously with the chiral 
phase transition [5]. However, at a nonzero chemical potential the results have so far been very limited, due to the 
severe problems lattice gauge theory encounters at nonzero /i. These difficulties arise from the fact that the action in 
the SU(3) gauge theory becomes complex- valued when fj, ^ [6], which in turn inhibits the importance sampling in 
the Monte-Carlo method. Although different computational schemes have been implemented [7-10] the problem with 
large-scale simulations at nonzero chemical remains unsettled. 

For N c = 2, the situation is much more favorable for numerical simulations. In this case, the fundamental rep- 
resentation is pseudoreal, and as a consequence the lattice measure remains real- valued even at fi ^ [11]. Since 
the fermion determinant is real, it is possible to perform lattice calculations for an even number of flavors using the 
standard Monte-Carlo algorithms. Although two-color QCD is worthwhile to study in its own right, the fact that it 
is a four-dimensional, confining non-Abelian theory makes it particularly interesting as a toy model for real QCD. 
Naturally, some caution is needed, since while there are some important commonalities between two- and three-color 
QCD, there are some striking differences. For example, hadron spectroscopy will be completely unrelated since for 
N c = 2 two quarks can form a color singlet. Likewise, the diquark condensate is a color singlet for N c = 2, while for 
N c = 3 it is not. Finally, in contrast to the color superconducting state of dense quark matter for N c — 3, diquark 
condensation for N c — 2 only gives rise to superfluidity. 
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Despite these differences, two-color QCD is still useful for testing general ideas and comparing analytical calculations 
with lattice simulations. Indeed, early lattice studies of two-color QCD yielded interesting results at both nonzero 
chemical potential and temperature [11]. In recent years, an extensive numerical program has been undertaken and 
a consistent picture of two-color QCD at nonzero T and \x has emerged [12,13]. For T = 0, as the chemical potential 
is increased, a colorless diquark condensate, (tptp), forms at [i c ~ m v /2, where m v is the pion mass. For \i > [i c , the 
condensate spontaneously breaks baryon number and gives rise to superfluidity This phase transition is most likely 
second order, at least for one and two flavors of staggered quarks for which the simulations have been performed. 
Alternatively, two-color QCD can be studied analytically using chiral perturbation theory. This effective theory 
describes the dynamics of the Goldstone modes associated with the spontaneous breaking of chiral symmetry and is 
valid for energies well below the chiral symmetry breaking scale. In particular, the critical baryonic chemical potential, 
H c = m n /2, is well within the effective range and indeed the onset of baryon density and diquark condensation at zero 
temperature is fully understood in this context at both the classical [14-19] and one-loop level [20]. The analytical 
predictions agree well with the lattice results, thus further strengthening our current understanding as well as showing 
the great utility of having both analytical and numerical tools available. 

Although the different phases and associated phase transitions at T = and fi ^ have been investigated in some 
detail, less is known about the phase structure of two-color QCD in the (/x, T)-plane. Analytical results employing 
renormalization group arguments for the chiral symmetry restoring phase transition at nonzero temperature but 
vanishing chemical potential were presented in Ref. [21], and the phase diagram at nonzero T and \x has been studied 
within a random matrix model [22]. On the lattice, simulations corresponding to Nf = 4 continuum flavors indicate 
that as T increases at a fixed moderate ratio x = 2^/mjr, the system undergoes a second order phase transition where 
(ipip) vanishes [12]. As this ratio increases, this line of second order phase transitions passes through a tricritical point 
and the phase transition becomes first order [12,23]. Recently the one loop calculation within chiral perturbation theory 
was extended to nonzero temperature [20,24], and the temperature dependence of the critical chemical potential for 
/i <~ m T /2 was calculated. Furthermore, the presence of a tricritical endpoint of this line of second order phase 
transitions was confirmed. The tricritical point was predicted to occur at T ~ [24] in qualitative agreement with 
the lattice results [12]. 

In this paper we further study the phase diagram at nonzero temperature and baryonic chemical potential by using 
a Landau theory combined with one-loop renormalization group arguments [25] . The Landau theory is fundamentally 
different from the chiral perturbation theory analysis in that it attempts to include all modes which become massless 
at the critical temperature. In addition to being a different approach, we are also considering a different range of the 
(/x, T)-plane than in Ref. [24]. Here, we study the chiral limit with a nonzero baryonic chemical potential. Since we 
consider the deep diquark phase a large temperature limit is of interest (in the chiral perturbation theory approach 
the low temperature limit was used for \i <~ m v /2). Our motivation is twofold. First of all, as mentioned above, the 
phase diagram of two-color QCD needs to be better understood, and one of our objectives is to provide a step in 
that direction. In particular, the emergence of the tricritical point at intermediate values of T and /i requires further 
explanation. Although the discussion will be concerned with the case of massless quarks, m q = 0, for simplicity, our 
results and predictions are expected to apply also for nonzero degenerate quark masses as long as Nf is even and 
/j 3> TOtt/2. Secondly, further comparisons between different analytical approaches and numerical simulations on the 
lattice will hopefully increase our understanding of both methods, and lay the groundwork for future attempts in full 
N c = 3 QCD at nonzero baryon density 

We assume throughout that the gauge group is SU (2) and that there are an even number Nf of massless flavors. In 
the next section, we discuss the pseudoreal representation of SU(2) and its implications for the N c = 2 QCD action 
and the symmetry breaking patterns. We then incorporate these properties into an effective Lagrangian. In Sec. Ill, 
we present our renormalization group analysis and discuss the order of the phase transition for various Nf. We first 
assume the absence of the chiral anomaly and later consider its implications. We end with our conclusions and an 
outlook in Sec. IV. 



II. TWO-COLOR QCD AND EFFECTIVE THEORIES 
A. The microscopic action 

The matter part of the QCD action with massless quarks and at zero chemical potential is 

Nf N, 

S = -iJ2 d ' x $f ( - i9A1(T a )l) r^jf = ~i £ / d'x , (1) 

f=i J f=i J 
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where the Minkowski metric is g^ u — diag(H ) and the color matrices are given by T a = a a /2, with a a the usual 

Pauli matrices. Defining = (1, (?) and a M = (1, — a), we have in the two-component formalism, 

Rewriting the action in Eq. (1) in this Weyl notation gives, after a partial integration, 



So = f d 4 x \ jf *» (d,8{ + igA;(T a ){) \) + (d,5j igA;(T a ){) & 

f=i J 1 



(3) 



Using the pseudo-reality relation in SU(2), a2<J a <^2 = —&a, where the superscript T denotes the transpose, the 
interaction term for A can be reorganized as 

\ jf {T a )l\) = A/TjA/ - -{\f*2)T a (*2\ f ) . (4) 

Since A/ transforms as 2 under a gauge transformation, cr 2 A/ transforms as £, i.e. like 2. Thus, by making the 
following flavor decomposition within the same color representation, 

^ f -\a 2 Xf f = N f + l,...,2N f ' ^ f -\\ f a 2 f = N f + 1, . . . , 2N f ' W 
we can write the action in Eq. (3) as, 

2JV> , 

So = ~i / d4x Qf - i9^T a ) a»Q f . (6) 

As can be seen from Eq. (6), the naive global symmetry S\J(Nf)i J xS\J(Nf)fix\J(l)B is in fact enlarged to SU(2Af) 
in two-color QCD, as a direct consequence of the pseudo-reality condition. 

However, this enlarged symmetry exists only as long as /x = 0. By introducing a common nonzero baryonic chemical 
potential /x, the action in Eq. (1) is augmented by an additional term, 

S ^S = S Q -nJ2 / d^frVi/ = So - n / d A x QBQ , (7) 
where B is the (2Nf)x(2Nf) diagonal baryonic charge matrix, 

Bs (o -°i) • W 

and and 1 are Nf x Nf matrices. The /x-dependent term is not an SU(2A/) singlet, since the block-diagonal form is 
only invariant under rotations that do not mix the different blocks. In other words, when fi ^ the global symmetry 
is reduced to the more familiar SU(A^)i,xSU(A^)ijxU(l)s. Furthermore, in the presence of a common quark mass, 
Tfiq , the symmetry is reduced even further, since the mass term explicitly breaks the axial symmetry. The global 
symmetry in this case is only SU(Nf)v xU(1)b- 



B. Spontaneous breaking of symmetries 

When n = m q = 0, we assume, as indicated by lattice simulations [12,13], that the chiral SU(2iV/) symmetry is 
spontaneously broken at low temperatures. A Lorentz and gauge invariant chiral condensate must be of the form [26] 

W oc QiE ij Qj , (9) 

with i and j the flavor indices and E a (2Nf)x(2Nf) antisymmetric matrix. Assuming a maximal flavor symmetry 
scenario [26], the standard form for E is given by 
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If the chiral condensate is nonvanishing, the global symmetry is spontaneously broken to the subgroup that leaves 
E invariant, namely Sp(2A/), and the symmetry breaking pattern SU(2Af )-^Sp(2A/) gives rise to (2N 2 — Nf — 1) 
Goldstone bosons. 

On the other hand, when fi > and the quark mass m q ^ 0, there is only a SU(A/)y xU(1)b symmetry. Although 
the chiral condensate does not break any of these symmetries, the lattice results show that the preferred ground-state 
configuration is still the chirally broken one in the region where [i < m^/2 [12,13]. However, at a baryonic chemical 
potential [i — m^/2, on-shell baryons are produced out of the vacuum, the chiral condensate begins to decrease and 
a diquark condensate appears: 



^ oc Q T £Q , (11) 



where £ is a (2A/)x(2A/) matrix given by 



I \ .„ T f -1 

I 



with I = I J q J ! (12) 



and I is an Nf x Nf matrix. A nonvanishing diquark condensate (m q ^ and fi > m n /2) breaks the SU(A^)y xU(l)s 
symmetry down to Sp(Nf) and gives rise to a Bose condensed phase. 

When the quark mass m q = 0, also vanishes. This situation requires some additional considerations for the 
diquark condensate. There is an interesting relationship, peculiar to two-color QCD, between the chiral and diquark 
condensates, at \i = 0. Consider the particular flavor transformation matrix Mi e SU(2A/), 

Mi = exp [(37ri/4)X] ; X=(°j j). (13) 

Under an SU(2A^) rotation, the chiral condensate transforms as 

E -» MiEMf = £ . (14) 

Hence, there is no distinction between the chiral and diquark condensate, and due to this identity we may as well 
define the symmetry breaking at /j = (and m q = 0) to be caused by a diquark condensate. There is no difference, 
and the breaking pattern remains SU(2A^)— >Sp(2A^). (For further discussion see [28].) 

C. The Landau theory at n — 

To model the m q — symmetry breaking, and its restoration as the temperature is increased for [i = 0, one 
can start from a Landau theory [21]. We assume that the Landau field can be parameterized by an antisymmetric, 
complex- valued (2Nf)x(2N f ) matrix $, $ — QQ T , that transforms under the SU(2A/) symmetry as $ — > M$M T 
with Me SU(2A/). The effective Lagrangian 1 which is renormalizable and invariant under the global SU(2A/) flavor 
symmetry is 

Lo = ^Tr[(d li &)(d"§)] -^-Tr[$t$] — Ai(Tr [$t$]) 2 - A 2 Tr[($ t $) 2 ] - c [Pf(*) + Pf(*+)] . (15) 
The last term in Eq. (15) is the so-called Pfaffian [27], 

Pf(<f>) = 2WN? E^ 1 )^- • • • > ( 16 ) 

' P 

where ik = 1, • • • , 2Af , and P denotes a summation over all permutations of {ii, . . . , i2N f } with (— l) p the sign of the 
permutation. This term ensures that the axial XJ(1)a is explicitly broken, since if c = the Lagrangian is actually 
invariant under the larger U(2A/) group. Stability of the classical potential requires that A2 > and Ai+A 2 /(2Af ) > 0, 
and by choosing m 2 < at T = the potential is minimized by a nonzero expectation value $ = Po^, with 
Pq = — m 2 / (8Nf\i + 4A2) for c = 0. This gives the correct symmetry breaking pattern SU(2A^-)^Sp(2A^), where we 
for definiteness took $0 to be proportional to £ instead of E. As discussed above, these are equivalent choices. 



1 We thank F. Sannino for access to his notes in which this effective Lagrangian was first written down. 
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D. The Landau theory at fi =fc 



We now introduce a nonzero baryonic chemical potential into the effective theory. The global symmetry SU(2iVf) 
is reduced to SU(A r /)ixSU(iV/)^xU(l)B and consequently there is a distinction between the chiral and the diquark 
condensate. Since we set = 0, the ground state has (tptp) ^ and (V>V) = at any fi > 0, breaking the symmetry 
spontaneously down to Sp(Nf)xSp(Nf). These aspects must of course be reflected in the effective Lagrangian when 
a chemical potential is added to it. A general discussion of Landau theories and the e-expansion at nonzero chemical 
potential is given in Ref . [29] . The important point to note here is that due to the pseudoreality of the representation, 
the order parameter is the same 2 as at [i = 0, namely QQ T . As in mean field theory, we assume that the quartic 
coupling constants remain unaffected by the introduction of /i, although we will later consider the possible influence 
of /i on the coupling c in Eq. (15). The dependence of the Lagrangian on /i is fixed uniquely by the symmetries and 
renormalizability. The transformation properties of the order parameter field is identical to that of the Goldstone 
field in chiral perturbation theory. Hence, the coupling of to $ enters through the covariant derivative [15] 

d v ^ d„$ - in5 v0 (B$ + $B) , 

d„&^ a„$t + i/iSv0 (B$t + $t B ) . (17) 

The Lagrangian in Eq. (15) then becomes 

L a ^ L = ^Tr [(9 M $ t )(0' 1 $)] + i^Tr [(a $)($ f B + B$ f ) - h.c] + yTr [(B$ + $B)(B$ f + $ T B)] 

2 

-^-Tr[$t$] _ Ax(Tr[$t$]) 2 _ A 2 Tr[($t$) 2 ] _ c [pf($)+Pf($t)] . ( 18 ) 

Given that $ transforms as <f> — > M$M T , the Lagrangian in Eq. (18) is only invariant under transformations that 
satisfy BM = MB, which reduces the global symmetry to S\J(Nf) L xS\](Nf) R x\](l) B for /i ^= 0. Moreover, from 
Eq. (18) it is straightforward to verify that the minimum of the classical potential is obtained when $ = 0o^, where 
(/)q = —(to 2 — i/j, 2 )/(8Nf\i + 4A 2 ) for c = and to 2 < 4^ 2 . This expectation value breaks the symmetry down to 
Sp(Nf)xSp(Nf), as required. It should also be noted that the choice <E>o = </>oE now corresponds to a local maximum, 
so that the effective model does distinguish between the two vacua. 



III. MELTING THE DIQUARK CONDENSATE 

As the temperature is increased, at a given /x > 0, we expect the diquark condensate to eventually melt and the 
symmetry to be restored. The object is to determine the order of this phase transition by studying the fixed points 
of the /3-functions in the e expansion. The effect of the chemical potential in the Landau theory is threefold: 1) it 
allows for terms linear in the time derivative, 2) it explicitly breaks the U(2iV/) symmetry, and 3) it induces a Bose 
diquark condensate. 

We will assume that T c at /i = is large, T 2 3> |to 2 |. Then as discussed in detail in [29] we can dimensionally 
reduce the Landau theory to an effective theory of the massless zeroth Matsubara modes. It is in this dimensionally 
reduced theory we will study the stability of the /3-functions. Because of the breaking of the U(2iV/) symmetry not 
all zeroth Matsubara modes become massless at T c when /j, > 0. At nonzero /j, the effective three dimensional theory 
describes the static zeroth Matsubara modes from the block diagonal modes in <!>. These are the massless modes at 
the phase transition. Explicitly, in the three dimensional theory the field <& is 




with Ai and A 2 Nf x Nf matrices given by, 



2 We consider values for the chemical potential which are not asymptotically large so that the effective degrees of freedom are 
given by bosonic diquark bound states and not by the quarks and gluons themselves. A more restrictive upper bound on /i 
is given by half of the mass of the lightest vector meson. At this value of the chemical potential, the formation of a vector 
condensate is expected [30-32]. 
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Ai = 



Pl2 





\—<f>lN f —4>2N f 



4>2N S 

o ) 



and A 2 = 



( <t>(N f + l)(N f +2) ■■■ <j>(N f +l)(2N f )\ 

-<f>(N,+l)(N,+2) ••• 4>( Nf+ 2)(2N f ) 



V -4>(N f +l)(2N f ) -(p(N f +2)(2N f ) 







(20) 



When the chemical potential is infinitesimally small, the masses of the zeroth Matsubara components of the off- 
diagonal blocks in $ are not large. In this particular situation, the phase transition is almost identical to the phase 
transition in the theory for /i = 0, at least if it is first order [21]. We therefore restrict /j, to be sufficiently large for 
this separation of mass scales to hold, although it should be emphasized that the matching to the nonzero T, n = 0, 
line in the phase diagram may be somewhat ambiguous. 

Our approach is to use renormalization group arguments in (4 — e) dimensions and study the existence of infra-red 
fixed points in the space of coupling constants for e — ► 1. Whether or not this extrapolation is valid is of course a 
difficult question to answer fully. However, this approach has been successful under many different circumstances, 
both in the condensed matter context (see e.g. [33,34] and references therein) and in the context of relativistic field 
theories [35,36]. 



A. Renormalization group flow in the anomaly-free case 

Having defined the model, we now analyze the renormalization group flow in the space of parameters Ai, A 2 , and 
c. As long as the quartic coupling constants are within their stability region, two situations can occur. Either the 
/3-functions have infrared stable fixed points, in which case there is a second order phase transition, or they do not, 
and the phase transition is then first order, induced by fluctuations. The obvious caveat here is that the coupling 
constants can be adjusted to lie outside the stability region, with the result that higher dimensional terms must be 
added to stabilize the theory, making the phase transition first order already at the mean field level. In addition, 
there could also be fixed points at the nonperturbative level. 

From the above discussion, the c — effective Lagrangian around criticality is 



L = ilr [(d v &)(ff'$j\ _ Ai (Tr [$t$]) 2 - A 2 Tr[($t$)' 



(21) 



where $ is defined as in Eq. (19). To study the existence of fixed points, we expand the Lagrangian in Eq. (21) 
to one-loop order around the expectation value <E> = an d calculate the j3- functions using the background field 
method. To more easily separate the one- loop corrections to Ai and A 2 , we will also label the constant 0o differently 
for different matrix elements: 



$0 -> i 



Pal 

fal 



(22) 



We then replace $ — > $0 + $ in Eq. (21) and split the fields into their real and imaginary parts, <j> mn — (<p$ n 
#ml)/\/2. Up to terms quadratic in the fluctuations, we find that 



, N f 



-Ai 



(d v ^f + (d^it) 2 + {dJ$ +Nf)(l+Nf) ? + (d^ +Nf){l+Nf) ) 

3N f /2 



kd = l 
Nf/2 

^k(N f /2+k) 



k<l 



k=l 



E *f 



k{Nj/2+k) 



k=N f +\ 



('-SI' 
E 
fe=i 



Nf/2 3N f /2 

k(N f /2+k) j\ ^k(N f /2+k) 
' ^k=N f + l 



N, 



N f 

-a, m x;{(A r) ) 



( ^ •)^2 , ^(^2 , ( Ar) ]2 , (Ai) 

VP(Nf+k)(N f +l)) + W(Nf+k)(N f +l) 



k,l=l 



2N f 



# 2 E {(^ ) ) 2 + (^ ) ) 2 > 



k < 1 k,l=N f + l 



k<l 
Nf/2 



k<l 



^{e{ 

k fe=l k 



J k(Nf/2+k) 



Nj/2 



k,l=l 



J k(Nf/2+l)Vl(Nf/2+k) f J k(Nf/2+l)ri(N f /2+k)~ r ^kl V( Nf /2+k)(Nj /2+J) ^kl Y(N f /2+k)(N f /2+1) 



J k(N s /2+k)l 



k<l. 
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3N f /2 

2 E 



,(r) 

k(N s /2+l) ( Pl(N f /2+fe) 



D k(Nf /2+l) ( Pl(N f /2+k) 



°kl C P(N f /2+k)(N f /2+l) 



5 /ci < P(N f /2+k)(N f /2+l) 



k,l=N f + l 



3N f /2 

E 

fc=JV / +l 



2 \2 
fc(Af//2+fe)-' 



where the classical potential Vq is given by, 



fe<i 



(23) 



(24) 



All the terms quadratic in (f> 2 and 4> 2 can be collected into a symmetric matrix M that can be diagonalized by 
an orthonormal transformation, M — > PMP T — Mjj. By integrating out the quadratic fluctuations and using 
dimensional regularization in = 4 — e dimensions, we find the effective action, 



V eS .(<Pa, <j> b ) = V - ^-^ (TrM 2 ) + 5MNf(<Pl + 4>lf + SXiNfiffc + <&) 



(25) 



where we have only kept the ultra-violet divergent terms and 8\\ and 5A2 denote the counterterms. Using the relation 
TrM 2 , = TrM 2 , it is straightforward to find the counterterms and hence the 0- functions. Using k as the arbitrary 
mass scale, we find 



01 = K-— 

OK 

02 = li- 



-eAi + 



-e\ 2 + 



2ir 2 



Xl{N f 2 - N f + 4) + 2X 1 \ 2 (N f - 1) 
5Aj (Nf - f) + f2AiA 2 



(26) 
(27) 



Consider now the different fixed points that can be inferred from Eqs. (26)-(27), i.e. the solutions to 0\ = 02 = 0. 
Since we are examining the long-range fluctuations, the fixed points A^ and X\ must be infra-red stable, implying that 
the eigenvalues of the stability matrix, 



, d0i/dXi 
~ { 002/ l dX 1 



80i/dX 2 
d0 2 /dX 2 



A2— 



(28) 



have to be positive. From Eqs. (26)-(27), we find the fixed points 



i : A* = 
iii : A* = 



A* =0 



^ 2 e 



Nf - N f + 4 ' 



ii: At = 0, X* = 



2^ 2 e 

m 1) 



5^/(^-1) -4 ' 



2ir 2 e(N f 2 - N f - 2) 
(N f -l)[5N f (N f -l)-4] ' 



(29) 



where the last one (iv) only exists for Nf > 4. In fact, when Nf = 2 it reduces to the third solution. From the 
eigenvalues of S, we find that the first and second cases are unstable for any number of flavors. The third case is 
unstable for > 4 but marginal when Nf — 2, i.e. one of the eigenvalues vanishes. The fourth solution, however, 
has a fixed point that is infra-red stable for Nf > 4. From this we conclude 

• For Nf = 2, the order of the phase transition is not determined by the renormalization group flow, and a higher 
order calculation is necessary to make a definite prediction. However, at /1 = 0, the temperature induced phase 
transition is first order for Nf — 2 and c = [21], so it is likely that this first order transition extends to small but 
nonzero values of \i. If the phase transition becomes second order as \i increases further, it should correspond 
to the 0(4) critical point (i.e. case iii above). 

• For Nf > 4, there is a stable fixed point. Thus, if this theory evolves from a point within the stability region 
X 2 > and Ai + X 2 /(2Nf) > 0, the phase transition in two-color QCD with more than four flavors, and without 
any influence from the axial anomaly, is likely to be second order. 



The renormalization group flow in the (Ai, A2)-plane for Nf = 4 is illustrated in figure f . 
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FIG. 1. The flow in the (Ai, A2)-planc for Nf = 4 and without the anomaly term. The classical potential is stable for A2 > 
and Ai + A2/(2iV/) > 0. The vectors shown have a length given by the magnitude of the eigenvalues of the stability matrix 
at the fixed point. The direction is away from the fixed point for negative eigenvalues and towards the fixed point for positive 
eigenvalues. The orientation of the vectors is given by the corresponding eigenvectors of the stability matrix. There is a stable 
fixed point and the phase transition is predicted to be of second order. For Nf > 4 the flow diagram has a similar topology. 
For Nf = 2 the upper right hand stable fixed point merges with the unstable fixed point on the Ai axis and becomes marginally 
stable. 



B. Adding the chiral anomaly 



The above conclusions must be reconsidered if the chiral anomaly is present around T c . Lattice results for N c = 3 
QCD indicate that although the effects of the chiral anomaly decrease, they do not vanish in the vicinity of the chiral 
symmetry restoring phase transition [37]. Since there is no reason to believe that the situation is drastically different 
in N c = 2 QCD, the anomaly can be expected to be present near T c and fj, = 0. At a nonzero chemical potential, 
there is an additional suppression of the anomaly [38] apart from the nonzero T effects. However, since our model is 
limited to the region \i < Aqcd, it is still likely that the effects of the chiral anomaly arc not negligible at the phase 
transition. 

For two flavors, the anomaly term is a mass- like operator and has the effect that only half of the number of the 
original real- valued fields in Ai and A2 become critical. At T = 0, there is an SU(2)xSU(2)xU(l) symmetry that 
is spontaneously broken to Sp(2)xSp(2), but since Sp(2)~SU(2), only the U(l)~0(2) symmetry is actually broken. 
Keeping only those fields in Eq. (18) which become critical at T c , the theory becomes 0(2)-symmetric and the diquark 
condensate breaks this symmetry completely. Thus, if the phase transition indeed is second order, it is characterized 
by 0(2) critical exponents. Since the phase transition at fi — can be second order as well [21], but with 0(6) critical 
exponents, there could thus be a change of the critical exponents as /x increases from zero. 

For Nf — 4, the anomaly term represents another quartic coupling. Repeating the calculation from the previous 
section with c ^ 0, we find the /3-functions, 



01 = k 

02 = K 

03 = K 



gAi 
dn 

d\ 2 
dn 
8c 
dn 



-eAi + 
-e\ 2 + 



16A? + 6A1A2 + 



04 



15A^ 



6AiA 2 + 



32 



-ec + ■ 



2c 



3Ai - A 2 



(30) 



When the fixed point c* vanishes, the solutions of course reduce to the results of the previous section. Surprisingly, 
there are no real-valued solutions to the three equations that satisfy c* ^ 0. Since there is then no fixed point, we 
conclude that the phase transition is first order, induced by fluctuations [39]. The presence of a first order phase 
transition for = 4 is also in agreement with the existing lattice results [12]. Note, however, (see Ref. [26,40]) that 
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the roles of the fundamental and the adjoint representations are reversed for the staggered fermions as used in Ref. 
[12]. Hence, this comparison is not fully justified. 

For Nf > 6, the critical behavior in the presence of the chiral anomaly is identical to the c = case. In renormal- 
ization group terms, the operator c is irrelevant, since its mass dimension in four Minkowski space-time dimensions is 
[c] = 4 — Nf, and the critical behavior is determined by Ai and A2 only. 

Finally, we stress again that in this approach we can not exclude first order phase transitions induced dynamically 
by higher order operators. Hence, even if there is a stable fixed point within the e-cxpansion, the transition may still 
be first order. 



IV. CONCLUSIONS 

In this work, we have studied the diquark phase transition in two-color QCD at zero quark mass and nonzero T 
and fj, using an effective Landau theory. To improve upon mean-field theory, we studied the renormalization group 
flow in 4 — e dimensions, and then extrapolated our results to e — 1. The conclusions to be drawn are the following. 
In the absence of the anomaly term and for Nf > 4, the phase transition could be second order, whereas the case 
of Nf = 2 is inconclusive at the one-loop level studied here. The effect of the anomaly on the phase transition has 
been shown to be substantial. In the presence of the anomaly, the Nf — 2 case exhibits a second order transition 
characterized by 0(2) critical exponents, which may change to 0(6) as fi — > 0. For Nf = 4, there is no infrared 
stable fixed points and hence the phase transition is first order. For Nf > 6, the phase transition may be of second 
order, since this situation is identical to the c = case. Of course, the statement for Nf > 6 is only true as long as 
N c = 2 QCD is asymptotically free. From the one-loop running coupling constant, this limits the number of flavors 

to Nf < iV/( m ax) = 11. 

Our analysis has focused solely on the case of massless quarks, and it is interesting to speculate on the consequences of 
adding nonzero masses. Defining x = 2fjb/m v oc fj,/ \/fn^, the phase transition discussed in this manuscript corresponds 
to the limit x — > 00. For a common, nonzero quark mass m q , the diquark condensate competes with the chiral 
condensate. Such a term alters the symmetry breaking pattern from S\J(Nf) L xS\J(Nf) R x U(l)s — >Sp(iVf)xSp(iVf) 
in the massless case to SU(A^)yxU(l)s — >Sp(Nf) when m q 7^ 0. For » m„ the first order phase transition for 
Nf = 4 and c ^ should only weaken, whereas the second order phase transitions for Nf > 6 presumably stay second 
order. When Nf = 2 the number of broken symmetries is independent of the introduction of a quark mass, since only 
the baryon number symmetry is broken. It is therefore possible that the phase transition is second order for all x. 
For Nf > 4, however, it is unlikely that the massless results should be representative at such a small ratio of /x and 
m n . Instead, the mass term should be incorporated explicitly into the effective Lagrangian, and the effects could very 
well change the order of the phase transition as x approaches 1. 

The presence of a tricritical point in the phase diagram of two color QCD has so far only been observed on the 
lattice for Nf = 4 flavors. Chiral perturbation theory [24] predicts that the tricritical point is present for any even 
number of flavors (Nf < Nf ( max ) = H)- hi contrast, the possibility of a fluctuation induced first order phase transition 
depends very sensitively on the number of flavors as well as the strength of the anomaly. By combining the results 
of future lattice simulations for different numbers of flavors with these complementary but quite different analytical 
approaches, the nature of the supcrfluid phase transition will certainly become more clear. 

In recent lattice investigations [41,42], it was found that the phase diagram of N c = 3 QCD at nonzero temperature 
and isospin chemical potential is very similar to that of iV c — 2 at nonzero temperature and baryonic chemical 
potential. In particular, the existence of a tricritical point was also demonstrated. Chiral perturbation theory also 
predicts the existence of a tricritical point for N c = 3 QCD at nonzero isospin chemical potential [24] . The approach 
followed in this work has been applied to this physically relevant case [29] and for Nf = 2 one marginal fixed point 
was found for the anomaly free case. 

By carefully combining the analytical and numerical results for a variety of parameters (including not the least of 
which is the case where Nf is odd), we should ultimately be able to grasp the complex phase diagram of N c — 2 QCD. 
Hopefully, this experience will also further enhance our knowledge of dense iV c = 3 QCD. 
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